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5 Field of the Invention 

The invention relates generally to the field of signal analysis, and more 
particularly, to a system and method for detecting the frequency, amplitude and/or phase 
of one or more tones comprised within an input signal. 

1 0 Description of the Related Art 

The discrete Fourier transform (DFT) is a popular tool for analyzing signals. 
However, before an input signal is transformed, it is quite often windowed with a 
windowing function. (It is noted that the action of capturing of a finite-length sequence 
of samples of the input signal automatically implies a rectangular windowing.) The 

15 transform Y of the wuidowed input signal will typically exhibit multiple scaled and 
shifted versions of transform function W, i.e., the transform of the window function. 
Each sinusoidal component of the input signal expresses itself as a pair of such shifted 
versions, one version shifted up to the frequency ^ of the sinusoidal component, and the 
other shifted down to frequency -f-. The positive frequency version is referred to herein 

20 as a positive frequency image, and the negative frequency version is referred to herein as 
a negative frequency image. When a sinusoidal component frequency is small 
compared to the sample rate, the positive frequency image and the negative frequency 
image for the sinusoidal component may overlap in frequency space. Similarly, when a 
sinusoidal component frequency is close to one-half the sample rate, the positive 

25 frequency image and the negative frequency image for the sinusoidal component may 
overlap. Furthermore, when two sinusoidal components have frequencies that are close 
together, their positive images and negative images may overlap. 



Atty.DktNo.: 5150-42800 



Page 1 



Conley, Rose & Tayon, P.C. 



Prior art techniques for tone estimation quite often focus on identifying the peaks 
in the magnitude spectrum |Y|. The peaks roughly determine the frequency of the 
corresponding tones. However, because of the cross-interaction of the images from other 
tones, or the negative frequency image from the same tone, the peak of a positive 
5 frequency image may be perturbed away from a purely scaled and frequency-shifted 
version of the template fimction W, Thus, parameter estimation techniques which 
compute parameters for a given tone based only on transform array values (i.e. DFT 
values) in the vicinity of a corresponding image peak may not produce accurate results. 
Therefore, there exists a substantial need for a system and method which could estimate 
10 tone parameters from the transform array with increased accuracy. 
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Summary of the Invention 

The present invention comprises various embodiments of a system and a method 
for estimating signal parameters (e.g. frequency, ampUtude and/or phase) of one or more 
sinusoidal tones present in an input signal. More particularly, one embodiment of the 
5 invention comprises a system and method for estimating parameters for a single tone 
based on a transform Y of the input signal. The input signal may be windowed with a 
window function w(n) and transformed into the frequency domain. The tone in the input 
signal expresses itself in the frequency domain as an additive combination of two spectra, 
one centered at the tone frequency and the other at the negative of the tone frequency. 

10 These two spectra are referred to herein as the positive frequency image and the negative 
frequency image respectively. The continuous-frequency transform W of the window 
fimction and the positive and negative frequency images have identically-shaped 
magnitude envelopes. Thus, a peak in the magnitude spectrum of the transform Y gives 
an initial estimate for the frequency and amplitude of the tone. Furthermore, the phase 

15 angle of the transform values in the neighborhood of the peak gives an estimate for the 
phase of the tone. The initial frequency, ampHtude and phase estimates may be used to 
compensate for the effect of a negative frequency image on the transform array in the 
frequency domain, especially in the neighborhood of the peak frequency. In other words, 
estimate values of the negative frequency image may be subtracted from the complex 

20 coefficients of the transform array in the neighborhood of the peak frequency. The 
resulting difference values may be used to compute improved estimates for the tone 
frequency, amplitude and phase. 

In one embodiment, a system may be configured to estimate signal parameters for 
one or more tones present in an input signal. The system may comprise an input for 

25 receiving an input signal, a memory, a processor and an output device, such as a display. 
The memory may store a software program which is executable by the processor. In 
response to execution of the software program, the processor is operable to perform the 
following operations. 
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(1) The processor may operate on samples of the input signal to generate a 
transform array. The transform array may include a positive frequency 
image and negative frequency image for each of the tones. The negative 
frequency images of the tones may distort or disturb observations of the 
positive frequency images of the tones of interest. Furthermore, one or 
more of the positive frequency images of the tones may also possibly 
disturb observations of the positive frequency images of the tones of 
interest. 

(2) The processor may identify locations of one or more magnitude peaks 
in the transform array. These frequency locations roughly locate the 
positive frequency images (and thus, the tone frequencies). 

(3) The processor may compute an initial frequency estimate, amplitude 
estimate and phase estimate for each of the one or more tones based on the 
transform array values in the neighborhood of a corresponding one of the 
peak frequency locations. These transform array values are complex 
numbers. The initial frequency estimate and amplitude estimate for each 
tone are determined based on the magnitudes of the transform array 
values. The initial phase estimate for each tone is determined based on at 
least one of the phase angles of the transform array values. These initial 
parameter estimates may incorporate errors since the positive frequency 
image of each tone is additively mixed with the other positive and 
negative frequency images. 

(4) The processor may correct the transform array values in the 
neighborhood of each peak frequency location by subtracting the effect of 
any interacting positive and/or negative frequency images due to other 
tones or the self-interaction due to the negative frequency image of the 
same tone. The positive and negative frequency images of the interacting 
(i.e. aliasing) tones are approximated using the initial frequency, 
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amplitude and phase estimates computed in (3) above. The difference 
values resulting from the transform corrections comprise a better 
approximation to the positive frequency image for each tone. 
(5) The processor may compute an improved frequency estimate, 
ampUtude estimate and phase estimate for each of the one or more tones 
based on the difference values, i.e. the corrected transform values, in the 
neighborhood of the corresponding peak frequency location. The 
improved frequency estimate and improved amphtude estimate for each 
tone are determined based on the magnitudes of the difference values in 
the neighborhood of the corresponding peak frequency location. The 
initial phase estimate for each tone is determined based on at least one of 
the phase angles of the difference values in the neighborhood of the 
corresponding peak frequency location. These initial parameter estimates 
may be more accurate than the initial estimate the effect of aliasing images 
on the image of interest have been compensate by the subtraction of step 
(4) above. 

The processor may transmit an indication of the improved set of signal parameters to the 
output device. For example, the improved set signal parameters may be displayed on a 
display device. Alternatively, the improved set of signal parameters may be forwarded to 
another system/device (or another software routine running on the same processor) for 
further processing. 

In step (1) above, the processor may window the input signal, and compute a 
discrete Fourier transform of the windowed input signal. The discrete Fourier transform 
may be implemented by a fast algorithm such as the FFT. 

The input signal may comprise one or more sinusoidal tones Xj, X2, x^ 
occurring at frequencies fj, f^, fL respectively. Each tone x^ expresses itself in the 
transform array as an additive combination of a positive frequency image of the form 
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(A,/2)exp(jejW(f-fO 
and a negative frequency image of the form 

(A,/2)exp(-jeOW(f+fi), 

where variable f denotes frequency, and W(f) is a continuous-frequency expression for 
5 the transform of the window ftinction w(n). Thus, the transform array comprises an 
additive combination of positive frequency images and negative frequency images 
corresponding to the one or more tones. Because the positive and negative frequency 
images may overlap with each other (especially when the tone frequencies are near zero, 
near one-half the sample rate, or near to each other), the frequency locations of magnitude 
10 peaks in the transform array may provide only a rough approximation to the tone 
frequencies f^. In other words, the observability of a given image may be adversely 
affected by the other positive and negative frequency images which overlap with the 
given image. 

The processor may identify the frequency locations of one or more magnitude 
15 peaks in the magnitude spectrum |Y(k)| of the transform array. In particular, the 
processor may search for magnitude peaks which exceed a magnitude threshold in a 
positive-frequency region of the transform array. A bin index value k^^^^ may be 
determined for each of the threshold-exceeding magnitude peaks. The bin index value 
k^x for each magnitude peak defines the bin index at which the corresponding magnitude 
20 peak is maximized. It is noted that the index k of the transform array is referred to herein 
as the bin index. 

The processor may compute a frequency estimate, an amplitude estimate and a 
phase estimate for each of the one or more tones based on a corresponding one of the 
magnitude peaks. The frequency estimate and ampUtude estimate for a given tone are 
25 determined from the magnitude values of the corresponding magnitude peak under the 
assumption that the magnitude peak is a shifted and scaled version of the window 
transform magnitude |W|. The center frequency of the magnitude peak determines the 
frequency estimate, and the size of the magnitude peak relative the window magnitude 
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|W| determines amplitude estimate. The phase estimate for a given tone is determined 
based on one or more the phase angles of the transform array coefficient (which are 
complex numbers) in the neighborhood of the corresponding magnitude peak. 

The frequency, amplitude and phase estimates for the one or more tones are used 
5 to estimate the positive and negative frequency images, and to subtract out the cross- 
interaction between images. More particularly, the processor may use the estimated tone 
parameters to correct the transform array values in the neighborhood of each peak 
frequency location. For a given tone, the processor may correct the transform array 
values around the corresponding peak frequency location by subtracting estimated values 
10 of any aliasing images. Aliasing images may include the positive and negative frequency 
images of tones other than the given tone, and the negative frequency image of the given 
tone. 

After correcting the transform array values, the processor may recomputed the 

tone frequencies, amphtudes and phases based on the corrected transform array values. 
15 Because the corrected transform array values more closely approximate the positive 

frequency images that the original transform array values, the recomputed parameter 

estimates may be more accurate. 

In one embodiment, the steps of correcting the transform array values and 

recomputing the parameter estimate may be performed repeatedly. When a termination 
20 criteria is achieved, the repetition may be terminated and final estimates for the signal 

parameters (e.g. tone frequencies, amphtudes and phases) may be transmitted to an output 

device (e.g. display screen). 

In one embodiment, the tone frequencies, amphtudes and/or phases may be used 

to decode analog and/or digital signal information contained within the signal. For 
25 example, the method may be used to more accurately identify the tones present in the 

input signal. Thus, the final estimates for tone frequencies, amplitudes and/or phases 

may be used to recover encoded analog and/or digital signals. 
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Brief Description of the Drawings 

A better understanding of the present invention can be obtained when the 
following detailed description of the preferred embodiment is considered in conjunction 
with the following drawings, in which: 
5 Figure lA illustrates a system configuration 100 for determining the signal 

parameters associated with one or more sinusoidal tones comprised within an input 
signal; 

Figure IB illustrates one embodiment for tone detection system 120; 
Figure 2 A illustrates one embodiment of tone detection system 120 comprising a 
10 computer-based measurement system, where signals generated by signal reception device 
SRD are presented to computer 102 through signal conditioning system 108 and data 
acquisition (DAQ) device 104; 

Figure 2B illustrates a second embodiment of tone detection system 120 
comprising a computer-based measurement system, where signals generated by signal 
15 reception device SRD are presented to computer system 102 through data acquisition 
(DAQ) device 104; 

Figures 3A-B presents a flowchart for one embodiment of a tone detection system 
according to the present invention; 

Figure 4 illustrates a windowing operation being performed on an input signal to 
20 generated a windowed input signal; 

Figure 5 illustrates the magnitude of transform array Y(k) for a typical windowed 
input signal comprising a single sinusoidal tone; 

Figure 6 illustrates a blowup of a generic magnitude peak 301 from the magnitude 
spectrum of Figure 5; 

25 Figure 7 illustrates the fact that the location and size of a magnitude peak may 

shift after correction, i.e. after subtracting cross-interaction terms due to ahasing images; 

Figure 8 illustrates the phase of transform array Y in the case where the window 
function is a Harming window; 
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Figures 9A-B illustrates one embodiment of a method for detecting signal 
parameters associated with one or more tones comprised with in an input signal; 

Figure 10 illustrates the magnitude of a transform array Y corresponding to a 
typical windowed input signal comprising three sinusoidal tones; 
5 Figure 11 illustrates three magnitude peaks Ul, U2 and U3 extracted isolated 

from the magnitude spectrum of Figure 10; 

Figure 12 illustrates the complex correction values D(k) for bin index values k in 
the neighborhood of a peak frequency location m^; 

Figure 13 illustrates an updated magnitude peak corresponding to original 
1 0 magnitude peak Uj. 



While the invention is susceptible to various modifications and alternative forms, 
specific embodiments thereof are shown by way of example in the drawings and will 
herein be described in detail It should be understood, however, that the drawings and 
15 detailed description thereto are not intended to limit the invention to the particular form 
disclosed, but on the contrary, the intention is to cover all modifications, equivalents and 
alternatives falling within the spirit and scope of the present invention as defined by the 
appended claims. 



20 
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Detailed Description of the Embodiments 



Figure lA 

Figure lA illustrates a system configuration 100 for performing signal processing 
5 on a signal comprising one or more tones. System configuration 100 may comprise a 
signal reception device SRD and a tone detection system 120. The SRD may coupled to 
receive a signal fi:om a device, unit under test (UUT) or a transmission medium 110, or 
any other system capable of transmitting a signal that may contain tones. The term 
"transmission medium" is used herein to refer generally to a device, unit under test 
10 (UUT) or a transmission medium 1 10 that may generate a signal including one or more 
tones. As used herein, the term "tone" includes a signal at a frequency, e.g., at a primary 
or single fi-equency, which may be contained within another signal. 

As shown in Figure lA, SRD may be coupled to a transmission medium 110. 
Transmission medium 1 10 may represent any of a variety of transmission media such as 
15 the atmosphere, free space, an optical fiber or fiber bundle, a communication bus (e.g. a 
network bus), a body of water or any other fluid, the earth, etc. In one embodiment, 
transmission medium 110 is the atmosphere, and signal reception device SRD comprises 
an antenna and a radio receiver. In a second embodiment, transmission medium 1 10 is a 
network bus connecting two or more computers, and signal reception device SRD is a 
20 network interface card/board. In a third embodiment, transmission medium 110 is an 
optical fiber, and signal reception device SRD comprises an optical sensor. As noted 
above, element 110 may be any of various devices or mediums for generating or 
transmitting a signal. 

Signal reception device SRD receives an input signal fi-om the transmission 
25 medium or device 1 10 and converts the input signal into a form suitable for presentation 
to tone detection system 120. The input signal may be electrical or non-electrical in 
nature. Signal reception device SRD may include analog-to-digital conversion hardware 
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to digitize the input signal. Alternatively, analog-to-digital conversion hardware may be 
comprised within tone detection system 120. 

In one embodiment, signal reception device SRD may comprise a measurement 
device such as a microphone, an accelerometer, a spatial displacement sensor, a strain 
gauge, a pressure sensor, a temperature sensor (e.g., a thermocouple), a radiation sensor, 
an optical sensor, etc, or any combination thereof. In another embodiment, signal 
reception device SRD may represent an array of transducers or measurement devices of 
one or more types. SRD may thus be any of various transducers or sensors for receiving 
a signal. 

Tone detection system 120 may couple to signal reception device SRD. Tone 
detection system 120 may be configured for detecting the fi-equency, amplitude and/or 
phase of one or more tones in the input signal. Tone detection system 120 may comprise 
a processor or central processing unit 140, memory 146, user input device(s) UID and a 
display device DD as shown in Figure IB. CPU 140 may be realized by any of a variety 
of computational devices such as a general purpose processor, a digital signal processor, a 
parallel processor, dedicated digital and/or analog circuitry, programmable gate array 
logic (e.g., an FPGA), etc., or any combination thereof Memory 146 may comprise any 
of a variety of memory devices such as random access memory (RAM) and/or read-only 
memory (ROM), as described further below. Tone detection system 120 may also 
include specialized data acquisition and/or signal conditioning hardware, interface 
hardware, etc., or any combination thereof 

Tone detection system 120 may comprise any of various devices, such as a 
programmable computer system, a computer-based system such as a VXI-based system, a 
PXI-based system, a GPIB-based system, a computer-based data acquisition system, or a 
dedicated test instrument, such as a dynamic signal analyzer, an oscilloscope or any other 
signal acquisition and/or analysis device. 

Tone detection system 120 may operate on samples of the input signal X 
generated by signal reception device SRD, and thus, may identify the frequency, phase 
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and/or amplitude of one or more tones in the input signal. The frequency, phase and/or 
amplitude of the one or more tones may be presented to a user through the display device 
DD or some other output device, and/or may be stored to memory for future use. 

User input device(s) UID may comprise a keyboard, a pointing device such as a 
5 mouse or trackball, a touch pad (such as those used in modem laptop computers for 
cursor control), a touch sensitive display screen, etc., or other input devices. Li one 
embodiment, user input device(s) UID may include use of a graphical control panel 
configured with various control icons such as buttons, knobs, sliders, switches, indicators, 
etc., or any combination thereof. A user provides input to tone detection system 120 
10 through user input device(s). Tone detection system 120 may manage a graphical user 
interface through display device DD and user input device(s) UID. 

Figures 2A and 2B 

Figure 2A and 2B illustrate exemplary embodiments of tone detection system 
15 120. As shown, tone detection system 120 may comprise a computer 102, a data 
acquisition (DAQ) device 104 coupled to the computer 102, and optionally a signal 
conditioning system 108 coupled to the DAQ device 104. Signal reception device SRD 
may comprise transducers, sensors, and/or receiving devices that couple to DAQ device 104 
through the signal conditioning circuitry 108. 
20 As shown, signal reception device SRD is configured and/or coupled to acquire 

signals from the transmission medium 1 10. The input signals acquired by signal reception 
device SRD may be optionally conditioned by the signal conditioning system 108 as shown 
in Figure 2A. The conditioned input signals may then be provided to DAQ device 104 as 
shown. Signal conditioning system 108 may connect to DAQ device 104 via one or more 
25 cables. 

Signal conditioning system 108 may comprise an external chassis 122 housing one 
or more signal conditioning modules 124 and optionally terminal blocks 126. Signal 
conditioning system 108 may be used to perform signal conditioning on field signals such 
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as the signals generated by signal reception device SRD. As used herein, the term "signal 
conditioning" may include one or more of amplifying, linearizing, limiting, isolating, 
filtering, switching and/or multiplexing field signals (e.g. transducer excitation), among 
other signal processing functions. Signal conditioning system 108 may advantageously 
5 reduce the level of noise in the signals transmitted to DAQ device 104. DAQ device 104 
may receive conditioned signals from signal conditioning system 108 as shown in Figure 
2 A. Altematively, DAQ device 104 may directly receive the input signal from signal 
reception device SRD as shown in Figure 2B. DAQ device 104 may operate to perform 
analog to digital (A/D) conversion and provides the resultant digital signals to computer 

10 1 02 for processing. 

Computer system 102 may include various standard components, including a 
processor or central processing unit (CPU) 140, system memory 146, non-volatile 
memory, one or more buses, and a power supply. DAQ device 104 may be a speciaUzed 
system for acquiring digital and/or analog signals from external devices. Thus, DAQ 

15 device 104 may include analog to digital (A/D) conversion circuitry and/or digital to 
analog (D/A) conversion circuitry. Examples of the DAQ device 104 include "E series" 
DAQ boards from National Instruments Corporation. DAQ device 104 may also 
comprise a computer-based instrument board, such as an oscilloscope, a digital multimeter 
(DMM), a dynamic signal analyzer, an arbitrary waveform generator, etc. 

20 In one embodiment, computer 102 may comprise input/output (I/O) slots into 

which DAQ device 104 may be coupled. In another embodiment, computer 102 may 
comprise a VXI (VME Extensions for Instrumentation) chassis and bus, a GPIB (General 
Purpose Interface Bus) interface card, a serial port or parallel port by which DAQ device 
104 may be coupled to the computer 102. 

25 Tone detection system 120, e.g., computer system 102, preferably includes at least 

one memory medium on which computer programs according to the present invention may 
be stored. The term "memory medium" is intended to include various types of memory or 
storage, including an installation medium, e.g., a CD-ROM, or floppy disks 104, a 
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computer system memory or random access memory such as DRAM, SRAM, EDO RAM, 
Rambus RAM, EPROM, EEPROM etc., or a non-volatile memory such as a magnetic 
media, e.g., a hard drive, or optical storage. The memory medium may comprise other 
types of memory as well, or combinations thereof. In addition, the memory medium may 
5 be located in a first computer in which the programs are executed, or may be located in a 
second different computer which connects to the first computer over a network. Li the latter 
instance, the second computer may provide the program instructions to the first computer 
for execution. Also, the computer system 102 may take various forms, including a personal 
computer system, mainfi-ame computer system, workstation, network apphance, Intemet 

10 appliance, personal digital assistant (PDA), television system, dedicated test or 
measurement instrument or other device. In general, the term "computer system" can be 
broadly defined to encompass any system having a processor which executes instructions 
fi-om a memory medium. 

The memory medium preferably stores a software program according to one 

15 embodiment of the present invention for detecting one or more tones in the input signal. 
More particularly, the software program may be operable to analyze the input signal to 
determine the fi-equency, phase and amplitude of one or more tones in the input signal. 

The software program may be implemented in any of various ways, including 
procedure-based techniques, component-based techniques, object-oriented techniques, or 

20 neural net based learning techniques, among others. For example, the software program 
may be implemented using ActiveX controls, C++ objects, Java objects, Microsoft 
Foundation Classes (MFC), or other technologies or methodologies, as desired. A 
processor, such as the host CPU, executing code and data fi-om the memory medium, or a 
programmable device configured according to a net list, may comprise embodiments of a 

25 means for determining the firequency, phase and amplitude of the one or more tones 
embedded in the input signal according to the methods described below. 

Various embodiments fiirther include receiving, storing, and/or transmitting 
instructions and/or data implemented according to the present invention upon a carrier 
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medium. Suitable carrier media include a memory medium as described above, as well as 
signals such as electrical, electromagnetic, or digital signals, conveyed via a 
communication medium such as networks and/or a wireless link. 

5 Figures 3 A&B - Aliasing Compensation Flowchart 

Figures 3A&B illustrate one embodiment of an aliasing compensation method for 
determining the frequency, amplitude and/or phase of a single tone present in the input 
signal. The method of Figures 3 A&B may be implemented by execution of a computer 
program stored on the memory medium as described above. 
10 In step 210, the CPU 140 may receive samples x(n) of the input signal provided 

by signal reception device SRD, and may multiply the input samples by a known window 
function w(n) to generate a windowed input signal y(n)=w(n)*x(n) as suggested by 
Figure 4. It is noted that the input signal samples may be received from a storage device 
(e.g. disk, CD-ROM) having been previously recorded/captured from signal reception 
15 device SRD. Alternatively, the input signal samples may be simulated samples generated 
by a simulator (e.g. a CPU executing simulation code). The present invention 
contemplates a wide variety of possible sources for the input signal samples x(n). 

The input signal is assumed to comprise a single sinusoidal tone in the presence of 
noise. Thus, the input signal may be modeled by the expression. 
20 x(n)=A*cos(G)on+0) 

-(A/2)exp(j6)exp(a)on)+(A/2)exp(-je)exp(-a)on), 
where 9 is the phase of the sinusoidal tone, A is the amplitude of the sinusoidal tone, 
(Oo=27rfo is the frequency of the sinusoidal tone, and n is a discrete time index. 

The window function w(n) may have any of a variety of forms. For example, the 
25 window function may be a rectangular window, a triangular window, a raised cosine 
window, a Hanning window, etc. 

In step 220, CPU 140 may perform a discrete Fourier transform (DFT) on the 
windowed input signal y(n) to generate a transform array Y(k), where k is a frequency bin 
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index which may range from 0 to N-1, or any interval of length N, where N is a positive 
integer. The transform array Y(k) may be modeled by the transform of the sinusoidal 
tone, i.e. 

Y(kHA/2)expae)W(f-fo)+(A/2)exp(-j0)W(f+fo), 
5 where W(f) represents the Fourier transform of the window w(n). It is noted that the 
relationship between frequency f and frequency bin number k is given by 

^fs*(k/N). 

where fg is the sample rate. The magnitude of the window transform W(f) typically has 
even symmetry and attains a maximimi at f=0. Thus, the function W(f-fo) attains a 

10 maximum magnitude at frequency f=fo, and the fimction W(f+fo) attains a maximum 
magnitude at frequency f^-fo- The first term in the expression above, i.e. 

P(f) = (A/2)expGe)W(f-fo), 
is referred to herein as the "positive-frequency image" since its center frequency occurs at 
the positive frequency f^. The second term in the expression above, i.e. 

15 N(f) = (A/2)exp(-j9)W(f+fo) 

is referred to herein as the "negative-frequency image" since its center frequency occurs 
at the negative frequency -f^. Thus, the transform array Y(k) includes a positive- 
frequency image and negative-frequency image which combine additively (in the sense of 
complex addition). The input signal may also include noise and/or other spurious tones. 

20 However, these are assumed to be insignificant for the embodiments described in 
connection with Figures 3A&B. 

If tone frequency fo stays away from zero or ^72, and/or, the sample size N is 
sufficiently large, the overlap between the positive and negative frequency images may be 
small, and thus, their individual identities may be apparent in the transform array Y(k). 

25 The magnitude ftmction |Y(k)| will thus exhibit two peaks which correspond to the 
positive and negative frequency images. The frequency locations of one of these peaks 
(i.e. the peak that occurs in the range of positive frequencies) may be used as an estimate 
for the tone frequency fp. 
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Conversely, if the tone frequency is close to zero or yi, and/or, the sample size N 
is sufficiently small, the positive-frequency image and negative frequency image may 
overlap significantly. Thus, their individual identities may not be apparent in the 
transform array Y(k). In other words, transform array Y(k) restricted to positive 
5 frequencies may be a poor approximation to the positive frequency image. Thus, the 
frequency location at which the magnitude ftmction |Y(k)| attains a maximum, when 
considered over positive frequencies, is only a crude initial approximation to the tone 
frequency 

Figure 5 is a plot of the magnitude of transform Y(k) corresponding to a typical 
10 windowed input signal y(n). Note that the transform Y(k) has a symmetry given by 
Y(k)-Y(k+N) for any integer k. Li particular, Y(-k)=Y(N-k). Thus, frequency bin 
numbers between N/2 and N may be mterpreted as negative frequencies. 

In step 230, CPU 140 may scan the DFT magnitude values |Y(k)| over the range 
of positive frequency bins to determine the bin index k which achieves the maximum 
15 magnitude. In other words, CPU 140 may select k^^^ as the integer bin index value k in 
the range from 0 to N/2 which maximizes the magnitude of Y(k). In addition, CPU 140 
may perform a comparison of |Y(k^^^-l)| and |Y(k^^+l)| to determine whether the second 
largest magnitude occurs at (k^ax-1) ^ Let denote the location of this second 

largest magnitude. Let a=|Y(k^J|, and let p= |Y(k2)|. 
20 It is noted that the maximum of magnitude fimction |Y(k)| considered as a 

ftmction of continuous frequency typically does not occur at the integer value k^^^, 
although it should occur somewhere in the interval between k^^^^^ and k2. Figure 6 
illustrates a blowup of the positive frequency magnitude peak 301 in the neighborhood of 
bin index value k^^. 

25 In step 240, CPU 140 may compute estimates and for the tone frequency fo 

and the tone amplitude A respectively based on the magnitude values |Y(k)| in the 
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neighborhood of the maximizing index k^^^ and an assumed functional forai for the 
window transform W(k). 

For example, in the case where the window function w(n) used in step 210 is a 
rectangular window, the window transform W(k) may be approximated by the expression 

W(k)=sin(7ik)/(7ik), Thus, the frequency estimate and real amplitude estimate may 
be computed according to the relations 

Ak=±p/(a4-p), 

TT Ak 

=a , 

sm{7r Ak) 

The plus solution for Ak is chosen if k2=kjjj^j^+l, and the minus solution for Ak is chosen if 

In the case where the window function w(n) used in step 210 is a Harming 
window, the window transform W may be approximated by the expression 

W(k)=sin(7rk)/[(7rk)*(l-k^)]. Accordingly, the frequency estimate and real amplitude 

estimate may be computed according to the relations 

Ak=±(2P-a)/(a+|3) 

sm{7r Ak) 
h=fsKIN 

Note that the plus solution for Ak may be chosen if k2=kj^ax"'"l5 the minus solution for 
Ak may be chosen if k2=kn^-l. 

A variety of window functions are contemplated. For some window functions 
w(n), it may be difficult to obtain a simple formula for the window transform W(k). In 
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these cases, values of the transform function W may be numerically approximated and 
used to compute the frequency and real amplitude estimates. 

In step 245, CPU 140 may compute an estimate 0^ for the tone phase using the 
phase angle of one or more of the complex values Y(k) in a neighborhood of k^^. In one 
5 embodiment, the phase of transform value Y(k^ J defines the phase estimate 4 , i.e. 

e,^angle{Y{k^J), 
where angle(z) denotes the principle angle of the complex number z. 

In a second embodiment of step 245, CPU 140 may interpolate the phase of Y(k) 
between and k2 to determine the phase estimate. For example, CPU 140 may 
10 perform a linear interpolation based on the phase of Y(k^J, the phase of Y(k2), and the 
value Ak. 

In other embodiments of step 245, CPU 140 may determine the phase estimate 0^ 
according to either of the expressions: 

^0 = angle{Y{floor(kQ))) or 

15 0,^angle{Y{ceiI(k,))), 

where floor(x) denotes rounding towards minus infinity, and ceil(x) denotes rounding 
towards plus infinity. 

As noted above, the transform array Y(k) is an additive combination of the 
positive frequency image and the negative frequency image, i.e. Y(k)=P(k)+N(k). 

20 Because the positive and negative frequency images may overlap (around DC and/or 
around Nyquist depending on the value of the tone frequency fo), the peaks appearing in 
the transform array Y(k) may be interpreted as disturbed versions of the corresponding 
images. However, given the estimates for tone frequency, amplitude and phase computed 
in steps 240 and 245, it is possible to compute the DC-aliasing and Nyquist-aliasing 

25 contributions of the negative frequency image on the transform array Y(k) in the 
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neighborhood of k^. By subtracting these aliasing contributions from the transform 
array Y(k), a better approximation to the positive frequency image may be obtained. 

In step 250, CPU 140 may use the phase estimate 9^ , the amplitude estimate , 

and the frequency estimate to compute the "DC-ahasing" contribution of the negative 
5 frequency image at frequency bins k in the neighborhood of k^j^. For example, CPU 140 
may compute estimated values A^^^ {k) of the negative frequency image according to the 
expression 

for bins k = floor{k^)-\-i-l, where i equals 0, 1, 2 and 3, and where floor(x) is the 
10 fimction which rounds x towards minus infinity. (It is noted this neighborhood of k^^^^ 
comprising four bins and starting at floor(ko)-l represents one of many possible choices.) 

In step 255, CPU 140 may use the phase estimate 9^ , the amplitude estimate , 

A. 

and the frequency estimate k^ to compute the "Nyquist-aliasing" contribution of the 
negative frequency image at the frequency bins k in the neighborhood of k^^^^. For 
15 example, CPU 140 may compute estimated values N^^^{k) of the negative frequency 
image according to the expression 

N,^^{k) = ^^xi^{-j9,)w[k-{N-kS) 

for bins k = floor{k^ ) + ^ - 1 , where i equals 0, 1, 2 and 3. 

In step 260, CPU 140 may compute estimated values P{k) for the positive 
20 frequency image according to the expression 

for the bin index values k in the neighborhood of k^^^. 
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attains a second- 



It is noted that the bin location kj^^^^ of the maximum magnitude for the function 
P(k) may not be the same as for transform array Y(k) as suggested by Figure 7. Thus, 
the parameter k^^ may be updated, i.e. set equal to the integer bin index k at which P(k) 

is maximized as indicated in step 265, and a may be set equal to P(k^) . Similarly^ 

5 parameter k2 may be set equal to the integer bin index k where P(k) 

highest value, and p may be set equal to Pik^) 

In step 270, CPU 140 may compute a second estimate k^^^ for the tone frequency 
and a second estimate A^^^ for the real tone amplitude based on the complex difference 

values P(k) generated in step 260. CPU 140 may use any of the methods described 
10 above in step 240 to determine these second estimates. Because the complex difference 
values P(k) more closely approximate the positive frequency image than the transform 
values Y(k) in the neighborhood of k^^, the second estimates may be more accurate than 
the first estimates. In other words, since the effects of the negative frequency image have 
been substantially reduced or removed, the new estimates computed in step 270 may be 
15 more accurate. 

In step 275, CPU 140 may compute an improved estimate 0^^^ for the tone phase 

based on the phase angle of one or more of the complex numbers P(k) in the 
neighborhood of the updated 1^^^,,. Any of the methods used to compute the phase 
estimate of step 245 may be used here to compute the improved phase estimate with the 

20 provision that P(k) substitutes for Y(k). 

In one embodiment, steps 250 through 275 may be iterated as many times as 
desired, or as many times as necessary to obtain convergence of the frequency, ampUtude 
and/or phase estimates. In each iteration of steps 250 and 255, the negative frequency 
image may be approximated in terms of the most recent estimates for the tone frequency, 
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amplitude and phase. For example, in a second iteration of step 250, the DC-aliasing 
contribution of the negative frequency image may be approximated by the expression 

A (2) 

Nik) = ~±^exp(~j§('^ )W(k + kf^ ) . 

After step 275, or after multiple iterations of step 250 through 275, CPU 140 may 
output the final frequency estimate, real ampHtude estimate and phase estimate to a user 
through display device DD or some other output device. Alternatively, these estimates 
may be stored in a memory for later use by some other signal processing device, or 
another software application running on CPU 140. 

The embodiments described above may generate estimates for the tone frequency, 
amplitude and/or phase even when the positive and negative images overlap significantly. 
For example, the tone frequency may be close to DC or one-half the sample rate, and/or, 
the size N of the DFT may be small. 



Banning Window 

In steps 250 and 255 described above, a phase estimate 9^ is used to compute 
respectively DC-aliasing and Nyquist-aHasing contributions of the negative frequency 
image to bins in the neighborhood of k^. In the Manning window embodiment, the 
phase estimate may be handled in different ways depending on whether aliasing 
compensation is being performed about DC or about Nyquist. Namely, for DC aliasing 
compensation, CPU 140 computes phase value according to the expression 

q>^=7t + angle{Y{k^)), 

where kjr= floor (k^) and =k^^^-\-Ak , and the DC aliasing contribution of the 
negative frequency image according to the expression 
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for bins k = floor(kQ) + i~l , where i equals 0, 1, 2 and 3, and where |x| denotes the 
absolute value of x. 

The form of the above expression for the phase estimate arises from the fact that 
the phase of Y(k) makes a jump of % radians between and k^^±l when the window 
function is a Manning window. Figure 8 illustrates this 180 degree phase jump in a plot 
of the phase of transform Y(k) for a typical sinusoidal tone which has been windowed 
with the Harming window. Note that the phase at the Nyquist frequency is not shifted 
with respect to the phase at k^^. Thus, for Nyquist-aliasing compensation, CPU 140 
computes the phase estimate (p^ according to the expression 

^0 -=angle(Y(k^)), 

i.e. without adding 180 degrees, and computes the Nyquist-aliasing contribution of the 
negative frequency image according to the expression 



2 

for bins k = floor {k^) + / - 1 , where i equals 0, 1, 2 and 3. See the source code appendix 
for a realization of the Banning window embodiment of the aUasing compensation 
method written in Lab View™. 



Detection of Multiple Tones 

In certain situations, the input signal may include multiple tones having different 
frequencies. Figures 9A&B illustrate one embodiment of a method for detecting the 
frequencies, amplitudes and/or phases of multiple tones in the input signal. It is noted 
that the method of Figures 9A&B may be implemented as one or more software programs 
stored in memory 146 and executable by CPU 140. 

In step 310, CPU 140 may receive an input signal x(n), and may apply a window 
w(n) to the input signal x(n) to generate a windowed input signal y(n)=x(n)*w(n). The 
input signal x(n) may originate from transmission medium 1 10, and may be presented to 
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tone detection system 120 through signal reception device SRD. However, the present 
invention contemplates a wide variety of source for the input signal samples x(n). For 
example, the input signal samples x(n) be may read from a memory medium (e.g. CD- 
ROM, magnetic disk, etc.) having been previously recorded/captured from transmission 
5 medium 110. Also, the input signal sample x(n) may be simulated samples generated by 
a simulator (i.e. a processor executing in response to simulation code). 

In step 320, CPU 140 may compute the DFT of the windowed input signal y(n) to 
obtain a transform array Y(k). 

L 

The input signal may be modeled by the expression x{n) = ^x-(n), where Xi(n) 

10 represents the i*^ tone of L tones in the input signal. The tone x^ is assumed to have the 
form 

Xi(n) = A^*cos((i>in4-9.) 
-(Ai/2)exp(]0i)exp(a)in)+(Ai/2)exp(-j9i)exp(-Q),n), 
where parameter 0)^=2^^ is the frequency of the tone x^, parameter is the real amplitude 
15 of the tone x^, and parameter is the phase of the tone x^. The input signal may also 
include noise and/or other spurious tones. 

The transform of the i*^ windowed tone yi(n)=Xi(n)*w(n) may be modeled as the 
sum of a positive frequency image 

P,(f)=(A/2)expOeOW(f-fi), 

20 and a negative frequency image 

N,(f)=(V2)exp(-j0i)W(f+fi), 
where W is a continuous-frequency expression corresponding to the transform of window 
w(n). (The positive frequency image has a magnitude envelope which is centered at tone 
frequency f^. The negative frequency image has an identically-shaped magnitude 

25 envelope which is centered at frequency -fj.) Thus, transform array Y(k) may be modeled 
by a summation of positive and negative frequency images 
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i=l 

f-fs*(k/N). 

If the tone frequencies maintain a sufficient mutual separation from one another, are 
sufficiently far from zero and f^/l, and the sample set size N is sufficiently large, the 
5 frequency support regions of the positive and negative frequency images may be 
essentially non-overlapping or minimally overlapping. Thus, each peak in the magnitude 
spectrum |Y(k)| may closely approximate one of the positive or negative frequency 
images, and the frequency location of the magnitude peak may accurately approximate 
the corresponding tone frequency f^. (Recall, the positive frequency images are centered 

10 on the tone frequencies). 

Conversely, if any of the tone frequencies get too close together, too close to zero 
or fs/2, or N is sufficiently small, the positive and negative frequency images may 
significantly overlap, and thus, a peak in the magnitude spectrum |Y(k)| may only poorly 
approximate its corresponding positive (or negative) frequency image, and the center 

15 frequency of the magnitude peak may be perturbed away from the corresponding tone 
frequency f^. Figure 10 illustrates the magnitude spectrum of a windowed input signal 
comprising three sinusoidal tones. It is assumed that each positive-frequency magnitude 
peak corresponds to one of the positive frequency images P^, and each negative- 
frequency magnitude peak corresponds to one of the negative frequency images Nj. 

20 In step 330, CPU 140 may scan the magnitude spectrum |Y(k)| to determine the 

frequency location of magnitude peaks occurring over the range of positive frequencies as 
suggested by Figure 11. In other words, CPU 140 may search for integer bin values m^ 
which correspond to local maxima of the magnitude spectrum when considered over 
integer bin values in the range from 0 to N/2. Let equal the maximal magnitude value 

25 for each peak, i.e. ai=|Y(mi)|. The local maxima may be subjected to a minimum 
magnitude test so that low-level noise peaks and signal side-lobes may be rejected. 
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In addition, CPU 140 may perform a comparison of the magnitudes |Y(m^+l)| and 
|Y(mi-l)| for each peak location m^ to determine whether the second largest magnitude for 
the corresponding magnitude peak occurs at k^^+l or k=mj-l . Let p^ denote the location 
of this second largest magnitude. Let represent this second largest magnitude value, 
5 i.e. PHYfe)|. 

In one embodiment, CPU 140 may identify positive-frequency magnitude peaks 
which satisfy a magnitude threshold relative to the largest magnitude peak. For example, 
CPU 140 may select positive frequency magnitude peaks that are more than X decibels 
below the largest positive-frequency magnitude peak, where X is a user selectable value. 
10 Figure 10 illustrates the magnitude peaks associated with three positive frequency images 
PI, P2 and P3 and the corresponding negative frequency images Nl, N2 and N3. 

In step 350, CPU 140 may compute for each tone x^, i=l, 2, 3, L, an estimate 

for the tone frequency f^ and an estimate A. for the tone ampUtude, These estimates 
may be computed based on the transform magnitude values |Y(k)| in a neighborhood of 
15 corresponding positive-frequency peak location m^, and an assumed functional form for 
the continuous-frequency spectrum W. 

In one embodiment, the window ftmction w(n) is a rectangular window. Thus, the 
continuous-frequency spectrum W may be assumed to have the form W(k)=sin(7ck)/(7rk). 
In this case, the frequency estimate f. and amplitude estimate A- for tone x^ may be 
20 computed according to the relations 

Aki=±p/(a,+p,), 

A. = — , 

sin(;7' A^.) 

k. = m- + Ak- , 
fi=fskJN. 
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The plus solution for Alq may be chosen if pi=mi+l, and the minus solution for Alq may 
be chosen if p^^m^-l . 

In a second embodiment, the window function w(n) is a Harming window. Thus, 
the continuous-frequency spectrum W may be assumed to have the form 

5 W(k)=sin(7rk)/[(7T;k)*(l-k )]. In this case, the frequency estimate f- and amplitude 
estimate A- may be computed according to the relations 

A,^a, ' (1-M/) 

sin(:;r M.) 

10 fi=fskfN 

The plus solution for Ak^ may be chosen if Pi==mi+1, and the minus solution for Ak^ may 
be chosen if Pj=mj-1 . 

A variety of window functions are contemplated. For some window functions 
w(n), it may be difficult to specify a simple formula for the spectrum W. In these cases, 
15 the values of W(k) may be numerically approximated and used to compute the Jfrequency 
and amphtude estimates. 

In step 355, CPU 140 may compute, for each tone x^, an estimate 9- of the tone 

phase 9j using the phase of one or more the transform array values Y(k) in the 
neighborhood of positive-frequency peak location m-. Any of the methods discussed 
20 above in the single tone embodiments may be used for the phase estimation of step 355. 

/V /V ^ 

Given the estimates for tone frequency , tone amplitude A- and tone phase 0. , 
the corresponding positive frequency image may be approximated by an expression 
such as 

A. 

p,{k)=^tMje,Wik-k,). 
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and the corresponding negative frequency image N; may be approximated by expressions 
such as 

NXk) = j-^^pi-A)Wik+k,) or 
iV, (k) = j- exp(-74 )Wik -(N-k, )) or 
5 N,(k)^j-exp(-jS,){w(k + k,) + W(k-(N-k,))}. 

In step 360, for each value of the index j running from 1 to L (i.e. the number 
tones), CPU 140 may compute the contributions of the other ahasing images on the 
transform array values Y(k) in the neighborhood of positive-frequency peak location nij. 
10 More specifically, for each value of the index CPU 140 may use the image 
approximations given above to compute a complex sum 

for bins k in the neighborhood of positive-frequency peak location m^-. In other words, 
the complex sum D(k) may include the estimated values at bin k of each positive 

15 frequency image other than Pj, and the estimated values at bin k of all negative frequency 
images. Figure 12 illustrates the complex values D(nij-1), D(mj) and D(mj+1) 
corresponding to a positive frequency magnitude peak Uj. 

In step 370, for each value of index j running from 1 to L, CPU 140 may subtract 
the sum D(k) from the corresponding DFT value Y(k) at each bin index value k in the 

20 neighborhood of positive-frequency peak location mj. The resulting difference values 
S(k)=Y(k)-D(k) comprise an improved approximation to the positive frequency image 
peak Pj, Figure 13 illustrates the magnitude of the difference values in the neighborhood 
of positive-frequency peak location m,. The subtraction in step 370 may operate to 
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reduce or remove the effects of the other positive and/or negative jSrequency images on 
the positive frequency image of the tone of interest. 

In step 375, CPU 140 may update the integer peak locations based on the 
magnitude of the difference values S(k). Because of the subtraction operation of step 
5 370, the magnitude peaks in the difference function S(k) may be shifted in frequency with 
respect to the corresponding peaks in transform Y(k). For each j in the range 1 to L, 
CPU 140 may examine the magnitude values |S(k)| in the neighborhood of peak location 
mj (i.e. the original peak location m^ computed above in step 330) to determine the integer 
bin index value of the new maximum magnitude. This bin index value becomes the 
10 updated value of peak location m^. The parameter may be updated as the new maximal 
magnitude, i.e. the magnitude of S(k) at new peak location m^. Similarly, CPU 140 may 
update the second-to-max peak locations pj and their corresponding magnitudes pj. 

In step 380, for each value of the index j running from 1 to L, CPU 140 may 

compute a second estimate //^ for the tone frequency fj- and a second estimate Af^ for 
15 the tone amplitude Aj based on the magnitudes of the complex difference values 

S(k)==Y(k)-D(k) in the neighborhood of updated peak location m^. CPU 140 may use the 

same (or similar) methods as those described above in step 350 to determine the second 

estimates. Because the complex difference S(k) values more closely approximate the 

positive frequency image peak Pj, these second estimates may be more accurate than the 
20 first estimates. In other words, since the effects of the other negative and/or positive 

frequency images have been substantially reduced or removed, the new estimates 

computed in step 380 may be more accurate. 

In step 385, for each value of index j running from 1 to L, CPU 140 may compute 

a second phase estimate 0^ for tone phase Gj based on the phase angle of one or more of 
25 the differences S(k) in the neighborhood of updated peak location nij. Any of the 

methods discussed above in the single tone embodiments may be used for the phase 

estimation here. 
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In one embodiment, steps 360 through 385 may be iterated as many times as 
desired, or as many times as necessary to obtain convergence of the frequency, amphtude 
and/or phase estimates. In each iteration of step 360, the positive and negative frequency 
images that contribute to the sums D(k) may approximated in terms of the most recent 
5 estimates for the tone frequencies, ampHtudes, and phases. 

After step 385, or after multiple iterations of step 360 through 385, CPU 140 may 
output final estimates for the real amphtude, phase and frequency of each tone Tj as 
indicated in step 390, These final estimates for the multiple tones may be presented to 
the user on display device DD or through some other output device(s). Alternatively, 
10 these estimates for the various tones may be stored in a memory for later use by some 
other signal processing device, or another software appHcation running on CPU 140. 

The embodiments described above may generate estimates for the tone 
frequencies, ampHtudes and/or phases even when the positive and negative frequency 
images of the tones overlap significantly. For example, the tone frequencies may be close 
15 to DC, close to one-half the sample rate, and/or close to each other. Overlap may also be 
due to spectral leakage when the size N of the DFT is small. 

Apphcations 

Embodiments of the present invention may be used in various applications. In 
20 general, embodiments of the present invention may be used in any system where it is 
desired to detect sinusoidal tones present in a signal, e.g., where it is desired to detect the 
precise frequency, amphtude and/or phase of the tones present in the signal. For 
example, an embodiment of the present invention may be used in a DTMF (Dual Tone 
Multi-Frequency) system for detecting tones present in a signal, such as a signal 
25 generated by a keypad of a telephone. Embodiments of the present invention are also 
contemplated for use in applications involving sonar, radar (e.g. Doppler radar), 
frequency-shift keying apphcations, mechanical systems analysis, etc. For example, the 
reflections generated by multiple moving objects in response to a radar pulse have 
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distinct frequencies dependent on their radial velocities with respect to the radar station. 
Thus, the frequencies of the reflections are usable for tracking the multiple moving 
objects. In another example, a mechanical system excited with a physical stimulus (e,g. 
an impulse) may manifest vibrations at one or more frequencies. The frequency, 
5 amphtude and/or phase of these vibrations may provide information to a system analyst 
about the nature of flaws in the mechanical system. Embodiments of the present 
invention may be used in a wide variety of applications, i.e. in any apphcation where it is 
desirable to identify one or more tones present in an input signal. The above-mentioned 
applications are merely representative examples. 

10 

Although the system and method of the present invention is described in 
connection with several embodiments, it is not intended to be limited to the specific 
forms set forth herein, but on the contrary, it is intended to cover such alternatives, 
modifications, and equivalents, as can be reasonably included within the spirit and scope 
15 of the invention as defined by the appended claims. 
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